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We extend tests of "Naik" and "fat link" improvements of the Kogut-Susskind 
quark action to full QCD simulations, and verify that the improvements previously 
demonstrated in the quenched approximation apply also to dynamical quark simula- 
tions. We extend the study of flavor symmetry improvement to the complete set of 
pions, and find that the nonlocal pions are significantly heavier than the local non- 
Goldstone pion. These results can be used to estimate the lattice spacing necessary 
for realistic simulations with this action. 



1 Introduction 



Because of the much larger computer time required to generate sample lattices with small 
lattice spacing in full QCD, improved actions are even more important to full QCD simula- 
tions than to quenched simulations. However, because of this same high cost of computing, 
the effects of improving the action have mostly been tested in the quenched approximation. 

For Wilson quarks, a Symanzik improved gauge action with a clover quark action Q 
and a renormalization group motivated improved gauge action with the conventional Wilson 
quark action]^ have been used in studies of full QCD thermodynamics. These studies showed 
significant changes in the physics relative to simulations with the conventional action. 

Studies of improvements with Kogut-Susskind quarks have been slower to start, possibly 
because it is less urgent to cure the o(a^) errors in this formulation than the o{a) errors in 
the conventional Wilson quark action. Kogut-Susskind quarks are preferred for studies of 
light quark thermodynamics because of their remnant exact chiral symmetry, but studies of 
quantities such as the equation of state have been severely limited by lattice effects[3. 



Two improvements to the Kogut-Susskind quark action have been investigated. Adding 
a third nearest neighbor term0 improves the quark dispersion relation and "free field" 
quantities such as the Stefan-Boltzmann law on the lattice. The effects of this improvement 
on the quenched hadron spectrum have been studied by the MILC collaboration 0, where it 
was found that it improves rotational symmetry as measured by the dispersion relation of the 
pions, but has only small effects on flavor symmetry breaking or spin splittings of the hadrons. 
Dynamical high temperature studies using the Naik action with four flavors of quarks have 
been done by the Bielefeld group |^. A second kind of improvement — averaging over paths in 
the parallel transport of the quark fields, or "fattening", does improve the flavor symmetry]^]. 
This can be understood as a suppression of the coupling to high momentum gluons which 
scatter quarks from one corner of the Brillouin zone to another [^], and this understanding 
has been used as a guide to further reducing flavor symmetry breaking by more extensive 
fattening |T0[. Recently, in|]lT| it was shown that in SU(2) quenched spectroscopy the local 
non-Goldstone pion becomes almost degenerate with the Goldstone pion after smoothing the 
gauge fields with RG "cycling". 



This work explores the effects of these improvements of the Kogut-Susskind quark action, 
together with a Symanzik improved gauge action, on simulations with dynamical quarks. 
One motivation for this is to make sure that the promising results from the quenched ap- 
proximation also apply to full QCD, and a second motivation is to begin developing a base 
of experience that can be used to estimate the quality of a physics simulation done with 
these actions at a particular lattice spacing, or alternatively to estimate the lattice spacing 
required for a simulation of some particular quality. We study improvement in rotational 
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invariance diagnosed through the pion dispersion relation, and flavor symmetry restoration 
from the zero momentum masses of the pions. In addition to the local non-Goldstone pion, 
the "712" , we look also at the nonlocal pions O, [1^ . No attempt is made at this stage to 



study whether these improved actions show better scaling of mass ratios. 



2 The simulations 



The gauge action used in most of these simulations was a Symanzik improved gauge action |]14, 
p!5| , p!^ , 0. This is the same gauge action used in the MILC studies of improvement for the 
quenched spectrum[||. 

Sg = (^] ( E (Pf^^) - 7^(1 + 0.4805a,) ^ (R,,) - ^0.03325 a, ^ (C,,^) 

\ -J / \x;fj,<u ZUUq ^.^^^ Uq x■,^l<v<a 

(1) 

where P is the standard plaquette in the /i, v -plane, and R and C denote the real part of 
the trace of the ordered product of SU(3) link matrices along 1x2 rectangles and 1x1x1 
paths, respectively. Here [3imp = lO/f?^, and a, = — 41og(Mo)/3.0684. Because these are 
exploratory calculations, no effort was made to tune uq. Instead, the value found in Ref. |Q 
at (3imp = 7.1, Uq = 0.8441 was used. 

The conventional one plaquette action was used in some simulations for comparison. 



Here P^onv = G/g"^- 

The quark action used the hopping matrix 2m + p with p parameterized as 



+ C3 {U^{x)U^{x + fi)U^{x + 2fl)6y^x+3fi — "backwards Naik term") 

+ E ( {UAx)U^{x + u)Ul{x + /i) + Ul{x - u)U^{x - u)U^{x - z> + ft)} Sy^^+f, 

— "backwards staple term" ) ) (3) 
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Here ci is the coefficient of the conventional single-hnk term, C3 is the coefficient of the third 
nearest neighbor (Naik) term, and is the coefficient of the "fat hnk", or staple term. 
Note that this parameterization of the fattening is designed to be convenient for full QCD 
updating, and differs from the parameterization used in Ref. p|. 

Four different valence quark actions were used in computing the spectrum. The conven- 
tional "one-link" quark action has Ci = 1 and 03 = ^3 = 0. The "Naik" action has Ci = 9/8, 
C3 = —1/24 and = 0. For the "fat-link" action we arbitrarily weight each staple with one 
half the weight of the single link path, and normalize so that the total weight of the one link 
path is one: ci = 1/4, C3 = 0, and W3 = 1/8. In Ref. it was found that the improvement 
in flavor symmetry is fairly insensitive to the exact weight of the staples. Finally, for the 
"one-f-fat-fNaik" action, we keep the relative weight of the staple to the one link term at 
one half, but normalize the total of the single link plus six staples to 9/8: Ci = (9/8)(l/4), 
C3 = —1/24 and W3 = (9/8) (1/8). Notice that we have not tadpole improved the quark 
action by including a factor of I/mq the third nearest neighbor coupling, as was done in 
Ref. 0. This would have been reasonable for the Naik action, but for the one+fat+Naik it 
would have made much less difference, since the bulk of the coupling to the nearest neighbor 
sites also comes from three link paths, which presumably are boosted by the same factor. 
Here the main effect would have been a renormalization of the mass. 



The gauge action and gauge force computations used a generalized gauge action code 



written by Tom DeGrand and Anna Hasenfratz ||18|| . For the inversions of the fermion ma- 
trix, the conjugate gradient algorithm was used. The products of links connecting to the 
third nearest neighbor sites and the sum of the single link plus staples were precomputed 
before starting the conjugate gradient, so the computer time per iteration was about twice 
that of the conventional one-link quark action. The fermion force can be computed by a 



straightforward generalization of the algorithm described in |]19|. The actions used here are 
all in the class of actions where the fermion matrix is 2m times the identity matrix plus a 
p operator which couples even sites to odd sites. For these actions the Hermitian positive 
definite operator M'^M is block diagonal, coupling even sites only to even sites. Thus, in the 
standard way, the pseudofermion action can be written without further doubling the number 
of flavors from four to eight as 

where 0e is a pseudofermion defined on the even lattice sites {x + y + z + t is even). (For 
simplicity, we sketch the force computation for four flavors of quarks with the hybrid Monte 



Carlo algorithm; adapting it to two flavors is done in the standard way||19||. Roughly 



' ■ ^ ^(MtM)-L-|0e) (6) 
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= (Xel^(MtM)lXe) (7) 

where 0e is restricted to even lattice sites and Xg is (M'^M)~^(/)e, also restricted to even sites. 
The diagonal (mass) term in M^'M doesn't contribute to the force, so by defining Xo on odd 
lattice sites to be p X^, this can be written as 

dM dM^ 

+ . (8) 

Except for a minus sign from the adjoint, the two terms can be treated identically. Now 
M depends on a particular link matrix U^{x) through every path in the quark action which 
contains this link. After projecting out the traceless anti-Hermitian part, the terms where 
the link is going backwards {W appears in M) just give a factor of two. So, with p X^ stored 
in the odd sites of X, the force computation involves 



• Compute Xg = (M^"M)~^0e on even sites. 

• Compute Xo = p Xe on odd sites 

• For each path in the quark action, for each link going in direction +fi, calculate a 
contribution to the "force" , or change in the momentum associated with the link: 

— Parallel transport X along the path from each end of the path to the site of the 
momentum being updated. This involves multiplying by U or W , and by the 
Kogut-Susskind phase factors ri^{x). 

— Take the traceless anti-Hermitian part of \F){B\, where F and B are X trans- 
ported from the forwards and backwards ends of the path. 

— Multiply by —1 if the link is an odd numbered link in the path, and by —1 if you 
are updating the momentum on an odd site. This takes care of the minus signs 
coming from the adjoint in the last line of Eq. ^. 

— Since the contribute an unwanted minus sign for every plaquette enclosed by 
a loop, multiply by —1 for each plaquette enclosed by a loop consisting of the 
path being considered and some standard path. For example, when the Kogut- 
Susskind phases are absorbed into the link matrices, the staple contribution to 
the fat link gets a minus sign relative to the single link connection. 

— Multiply by the coefficient of the path, which includes a minus sign for paths that 
are part of the "backwards" part of p. 



Some optimizations are possible, since several different paths passing through a site may 
share the same forwards or backwards parts and vectors transported to one site and used 
in the force computation there can be subsequently transported to the next site along the 



4 



path. The timings of the various parts of the computation depend on many factors such 
as the machine used, the lattice sizes, and how carefully each part is optimized. As an 
example, on an 8^ lattice on an alpha workstation with this improved action the gauge 
force computation requires 24.6 seconds, the fermion force computation 5.4 seconds and the 
conjugate gradient 0.183 seconds per iteration. About 400 conjugate gradient iterations were 
required at /3imp — 7.5 and arriq — 0.015. 

Since getting good acceptance in the hybrid Monte Carlo algorithm is a good test of the 
program, some four flavor results on small lattices were done first. While these were intended 
just to be test runs, because they had large lattice spacing and larger quark mass than the 
two flavor runs they ended up showing the largest effects on the pion dispersion relation, so 
some of the four flavor results for the pion dispersion relation are reported. Two flavor runs 
with the improved action were done at f3imp = 7.1, 7.3 and 7.5. Some of the conventional 
action results needed for comparison can be found in the literature, but some conventional 
action lattices were generated at Pconv = 5.5. This was done to get a reasonable match of 
lattice spacing and lattice size to the Pi^p — 7.3 improved action runs and to get the masses 
of all of the pions. For all of the sample sets that were generated, the meson spectrum was 
calculated with both the conventional and improved valence actions. In some cases valence 
actions where the Naik and the fattening were turned on separately were also used. Several 
local meson propagators were calculated at nonzero momenta, but the Goldstone pion gave 
by far the most accurate results, and our analysis of rotational symmetry will concentrate 
on these propagators. 



3 Results 



The largest data sets are two flavor runs at Pimp = 7.1, 7.3 and 7.5, together with conventional 
action runs at Pconv — 5.5 chosen to match the Pimp = 7.3 runs. At Pimp = 7.1 a single mass 
value was run, at rriT^/mp = 0.51 for the conventional action and niT^/mp = 0.57 for the 
fat+Naik action. At the two larger Pimp values, two quark masses were used which bracket 
niTr/mp = 0.55, so it is possible to interpolate to this arbitrary value for comparing the 
different actions. Also, to compare the lattice spacings of the different runs we can use the 
P2 mass interpolated to rriT^/mp = 0.55 as an arbitrary length standard, and one which is 
more accessible than an extrapolation to the physical light quark mass. (The p2 is used 
instead of the p since at large lattice spacings it often gives a better signal.) 

Table I contains the pion masses, including both local and nonlocal pions at zero mo- 
mentum, and the Goldstone pion at momenta ^(0,0,0), ^(0,0,1) and ^(0,1,1). The p2 
mass is also tabulated, as it will be used to set the scale. 

The absence of rotational symmetry on the lattice leads to meson dispersion relations that 
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differ from the E{kY 




expected in the continuum. The dimensionless parameter 



"speed of hght squared" , which should be one, is used to parameterize the deviations from 
the correct dispersion relation, 



Table II shows for the Goldstone pion at the two smallest lattice momenta in selected 
runs. While the two flavor runs with smaller quark masses are relevant for estimating the 
lattice spacing required for physical simulations, improvement of the dispersion relation is 
naturally most dramatic in the runs where it is worst with the unimproved action, namely 
at large lattice spacing and large quark mass. Therefore Table II also contains a four flavor 
run at quark mass amq = 0.1, corresponding to m^/mp = 0.76 or 0.82 with the conventional 
and improved actions respectively. Ideally, all of these numbers would be quoted at the 
same physical momentum. A crude approximation to this was achieved by increasing the 
spatial lattice size in units of a as /5 was increased, hoping to keep the physical lattice size 
constant. The extent to which this was achieved is shown by expressing the unit of lattice 
momentum, 2n/L in units of the p2 mass interpolated to m-^/mp = 0.55, and is included in 
this table. Since involves the difference of the pion energy at zero and nonzero momentum, 
a jackknife analysis was done to estimate the statistical error on at (3 = 7.3 and 7.5. Errors 
estimated from the jackknife analysis were not significantly different than errors estimated 
naively from the statistical error in the nonzero momentum pion energies. The improved 
action = 7.3 and conventional action G/g"^ = 5.5 results in Tables I and II contain the 

non- Goldstone pion masses and speed of light for the Goldstone pion with the fat-link and 
Naik improvements turned on separately. Within the fairly poor statistics, this is consistent 
with the quenched result of Ref. that, as expected, it is the Naik term that is responsible 
for the improvement in the dispersion relation. 
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Table 1: Summary of vr and p2 masses from selected runs. The top two lines of the 
table show the format for tabulating the masses. The first column is the gauge coupling, 
where the subscript "C" indicates the conventional one-plaquette action and the subscript 
"I" indicates the Symanzik improved gauge action. The second column is the quark mass. 
The third column is the spatial lattice size, and the fourth column is the valence quark 
action. In this column "O" stands for the conventional "one link" action, "OFN" for the 
improved "one-link + fat + Naik" action. "OF" and "ON" stand for the "one-link + fat" 
and "one-link -|- Naik" actions respectively. The fifth column is the p2 mass. The remaining 
columns are pion masses. The first four columns of pion masses correspond to distance 0, 1, 
2 and 3 operators respectively, all at zero momentum. There are two pion operators at each 
distance, shown one above the other. These operators have the "spin ® flavor" structures 
shown in the table header [Q. The Goldstone pion is the upper pion in the distance zero 
column. Goldstone pion energies at nonzero momenta are in the last column, labelled by 
the momentum in units of 2tt/Ns. 
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Table 2: "Speed of light" for the Goldstone pion. The first three columns are the gauge 
coupling, quark mass and valence quark action. The fourth column is the unit of lattice 
momentum in units of the arbitrary mass scale nip at mT,/mp — 0.55. The last two columns 
are obtained from the lowest two nonzero lattice momenta. 
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Another symmetry which is broken on the lattice is the flavor symmetry. The pion masses 
with Pimp = 7.3 or Pconv = 5.5 with various valence quark actions in Table I show that the 
fattening of the link improves the flavor symmetry while the Naik term has little effect, in 
agreement with earlier studies in the quenched approximation!^ HOI- Here we also see that 
the improvement extends to the one-link nonlocal pions. We can parameterize the splitting 
of the pions by the dimensionless quantity 

Here is one of the non-Goldstone pion masses, Mq is the Goldstone pion mass, and Mp is 
one of the local p masses. Since the p masses are nearly degenerate, it makes little difference 
which one we use. In our analysis we used the p2 mass because it is usually estimated more 
accurately. This measure of flavor symmetry breaking differs from the quantity used in [||. 
Here we use the squared masses since they are approximately linear in the quark mass and 
we will do some interpolating and extrapolating. Also, empirically this S2 is fairly insensitive 
to the quark mass. The squared Goldstone pion mass is subtracted in the denominator to 
give more sensible behavior at large quark mass. At large amq all of the meson masses are 
becoming degenerate, but it is still possible to use the difference between the Goldstone pion 
mass and the rho mass as a natural scale for flavor symmetry breaking. The local non- 
Goldstone pion is the easiest to measure, and one for which results with the conventional 
action are available in the literature. Figures |I] and ^ show 62 for some of the pions as a 
function of lattice spacing, both for this improved action (both the dynamical and valence 
quark actions are improved), and for the conventional action. We show 62 using the local 
non-Goldstone pion (7075^7075), which is consistently the lightest non-Goldstone pion, and 
the distance three 7075 1 pion, which is usually the heaviest. Results for the other pion 
operators can be filled in using the values in Table I. Here we have plotted results for both 
conventional and improved dynamical actions in the same plot, with the valence actions in 
separate plots. Thus, to compare a completely conventional simulation (both valence and 
dynamical actions) with a simulation with both valence and dynamical actions improved, 
one should compare the burst and diamond in Fig. |1] with the octagons and squares in Fig. ^ 
In these figures all of the data have been interpolated to mT^/nip = 0.55 except for the points 
at = 7.1. 

Clearly, the fat-|-Naik action shows a significant improvement of the flavor symmetry 
breaking over the conventional action. This improvement is due to the "fattening" and not 
due to the Naik term which only affects the dispersion relation. The nonlocal pions have 
significantly larger splittings than those of the local H2. Although their splittings get reduced 
along with the 112 splitting, it is evident that the reduction is smaller. However, as expected, 
all the pions are approaching the same mass as the lattice spacing gets small. From our 
results, it is clear that one has to consider the nonlocal pions before making any judgment 
about the quality of the improved action. 
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Figure 1: 82 for the conventional (One Link) valence action. Points with improved dynamical 
action are at /Sj^p = 7.1, 7.3 and 7.5. All the points are interpolated to myr/mp = 0.55 except 
for the ^imy — 7.1 points for which mT^/nip = 0.51. The octagons correspond to the local 
non-Goldstone pion which has the smallest 82 and the squares correspond to the three-link 
7o75 ® 1 pion which has the largest 82- Results with the conventional dynamical action at 
f^conv = 5.5 are also shown, where the burst is from the local non-Goldstone pion and the 
diamond from the three-link pion. The diagonal lines are crude extrapolations to smaller 
lattice spacing, where the dotted hne is the one-|-fat-|-Naik valence action local pion, the 
dot-dashed line the one-|-fat-|-Naik valence action three link pion, the long dashed line the 
conventional valence action local pion, and the short dashed line the conventional valence 
action three link pion. 
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Figure 2: 82 for the fat+Naik valence action. The dynamical actions and the symbols used 
in plotting are the same as in Fig. |l]. All the points are interpolated to m^/mp = 0.55 except 
from the the (3imp = 7.1 points for which mT^jnip = 0.57. 
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Furthermore, it is evident that part of the flavor symmetry restoration is due to the 
improvement of the gauge action. This is also supported by the quenched spectroscopy 
results presented in[^]. In the case of the one link action (Fig. |l|), the improvement of the 
gauge action significantly reduces both the local and the nonlocal pion splittings. The local 
pion splittings are affected more than those of the nonlocal. In the case of the fat+Naik 
action (Fig. |^), we also see a reduction of all the splittings. However, the local pion splittings 
are improved very little while the nonlocal get a major improvement. From the above 
observations one may conclude that the "fattening" of the fermion action is less effective 
for the the nonlocal pions which they get improved mostly because of the the gauge action 
improvement. Further "fattening" may be needed in order to significantly improve the 
nonlocal pions. 



4 Conclusions 



The first conclusion is that, as expected, the improvements in rotational symmetry and flavor 
symmetry seen in the quenched spectrum also appear in full QCD spectrum calculations. 
The most serious problem for accurate simulations of full QCD with Kogut-Susskind fermions 
is the flavor symmetry breaking. While breaking of rotational symmetry is also serious at 
large lattice spacings, the light meson dispersion relations are in much better shape than the 
flavor symmetry. In particular, in our simulations with jSimp = 7.5 the speed of light for the 
pions is quite acceptable, while the flavor symmetry is still seriously broken. 

We have extended the earlier studies of flavor symmetry breaking to include the nonlocal 
pions as well as the local non-Goldstone pion. The splittings between the nonlocal pions 
and the Goldstone pion are also reduced by improving the action, but they are significantly 
more massive than the local non-Goldstone pion. Worse, the amount of flavor symmetry 
improvement from the fat link used here is less for the nonlocal pions than the local one. 
This is seen in Figs. |T] and 0, where we have drawn some lines crudely indicating S2 for the 
local pion and the 3-link 7075 ® 1 pion, which is consistently the heaviest pion. Note that 
the slope of the line for the local pion was reduced by a larger fraction than was the slope 
for the 3-link pion when the valence action was improved, i.e. in going from Fig. |I| to Fig. ^ 
Unfortunately, in a realistic simulation we would want all of the pions to be well behaved, so 
in estimating the lattice spacing required for the desired quality simulation we must consider 
the nonlocal pions. 



An important use of action tests such as this work is to estimate the lattice parameters 
required to investigate real problems in QCD. An important example, and the example which 
partially motivated this work, is the study of the high temperature QCD phase transition. 
Interesting questions here include an accurate determination of the equation of state [^], and 
the question of whether the addition of a strange quark changes the nature of the crossover 
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or phase transition [^]. To be sure that we are doing an accurate simulation in the crossover 
temperature range, we would like to know that our simulation is under control on both the 
high and low temperature sides of the crossover. The addition of the Naik term does an 
excellent job of improving the quark dispersion relation, and hence we expect to improve the 
equation of state in the high temperature phase^]. However, the harder part of the job is 
to get the physics right on the low temperature side. Roughly, there the physics is a gas of 
thermally excited pions. To reproduce this physics, we need the pion dispersion relation to 
be good up to energies large compared to the temperature, and the pion masses to be nearly 
degenerate — that is, for flavor symmetry to be restored. For concreteness, suppose that we 
are interested in the effects of adding a strange quark. Viewed from the low temperature side 
of the transition, this adds kaons with mass 500 MeV to the physics. To isolate this effect, 
we should require that all of the pions have masses small compared to this. For example, 
we might demand that when the Goldstone pion has its physical mass of 140 MeV, that 
the non-Goldstone pions have masses less than one half the kaon mass, or 250 MeV. This 
requires 62 < 0.075. The for the various pions are fairly insensitive to quark mass. ( 
In the 10/ gj = 7.5 runs they are equal within errors at arriq = 0.015 and arrig = 0.030 for 
both the one-link and one-fat-Naik valence actions.) Neglecting dependence of 62 on quark 
mass, it is straightforward to find the necessary value of arup from the lines in Figs. |l] and H, 
and from there to figure out the number of time slices needed to simulate a temperature of 
150 MeV. Requiring that the heaviest pion mass be less than 250 MeV when the Goldstone 
pion is at 140 MeV will require about 20 time slices with the conventional action, and about 
16 with the one+fat+Naik action. (Had we used the local pion in this exercise, the results 
would have been 16 and 10 time slices respectively.) Sixteen time slices is a disappointingly 
small lattice spacing, so it is important to extend these tests to actions which might further 
reduce flavor symmetry breaking, such as the action discussed in Ref. |T^. 
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